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1. INTRODUCTION 

Currently, the parallel transmission networks are widely used in the electrical power systems. 
Therefore, a fast and reliable protection is very much needed in such aspects as rapid fault detection and 
accurate estimation of the location fault will reduce errors, and assist in the maintenance and restoration 
services to improve the continuity and reliability of electric supply. Therefore, parallel transmission lines 
require more special consideration in comparison with a single transmission line, due to the effect of mutual 
coupling on the parallel transmission line. It also must conform to the standard IEEE.STD.114 2004 [1]. One 
main advantage of parallel transmission is the availability of the transmission lines during and after the fault. 

Fault location problems in the parallel transmission lines have been widely researched. Many 
diagnostic approaches have been proposed in the literature, including fault location based on the amount of 
electricity, which includes one-terminal method [2], [3], two-terminal method [4], [5], traveling wave 
analysis method [6], [7], resistance measurement method [8], and the determination of fault location 
estimation using wavelet transform [9]. 

This paper proposed discrete wavelet transformation using the Clarke’s transformation to determine 
the fault location estimation and classification on the parallel transmission lines. This study presents a 
different approach, which is based on the Clarke’s transformation known as alpha-beta transformation, which 
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is a transformation of a three-phase system into a two-phase system [10], [11], where after the result, the 
Clarke’s transformation is then transformed into discrete wavelet transform. 

Wavelet transform is an effective tool in analyzing the transient current and signal associated with 
faults, both in frequency and in the time-domain [12], [13] and is ideal for dealing with non-stationary signal. 
This can improve the accuracy, reliability of the detection and classification of power quality disturbance 
[14], and features can be applied to determine the fault location estimation [15]. 

The proposed approach combines the decomposition of electromagnetic wave propagation modes, 
using the Clarke’s transformation of signal processing, given by the discrete wavelet transformation based on 
the maximum signal amplitude (WTC) ° to determine the intrusion time. This work made extensive use of the 
simulation software PSCAD/EMTDC [16] which resulted in the fault of the simulation of the transient signal 
transmission line parallel to the number of data points (10°). For one kind of fault, these data were then 
transferred to MATLAB with the help of Clarke’s transformation to convert the three-phase signal into alpha 
and beta signals. The signals were then transformed into several mother wavelets [17] such as Db4, Sym4, 
Coil4 and Db8 which were manipulated for comparison in terms of time and the distance estimation fault in 
the parallel transmission line. 


2. BASIC PRINCIPLE OF CLARKE’S AND WAVELET TRANSFORMATION 
2.1. Clarke’s Transformation 

Clarke's Transformation (aß0) is a useful analytical approach to complement symmetry components 
(0, 1, 2). Clarke’s transformation can overcome some symmetry component drawbacks such as the 
calculation of the circuit transient phenomena. Clarke's Transformation (aß0) is the transformation method 
that contains the elements of a 3 x 3 matrix, containing matrix element in the form of real, whereas the 
symmetry component contains matrix components in the form of real and complex numbers. A three-phase 
current which has a digital representation is assumed to have the form [18], [19]. Therefore, the above 
components can be formed into a matrix [20] 


1 1 
E 2225 se 
a « ig (n) 
i r a 2 P 
İmode=İapo=Clapc =; 0 = -2j|x ip (n) (1) 
22 get 
2 2 


where C is the well-known transformation introduced by Edith Clarke [21]. 
2.2. Fault Characterization in Clarke’s Transformation 
2.2.1. Single Line of Ground Fault Phase A to G 
A suppose for a line to ground fault (AG), assuming the grounding resistance is zero, then the 
instantaneous boundary conditions will be: 


Ip = Íe =OandV, =0 (2) 


Then, the boundary condition instantaneous will be: 


la —2/3 I, 

Ip 0 

Ip} = | -1/3 Ig (3) 
L -2/3 I, 


2.2.2. Line to Line Fault in Phase A-B 
A suppose for the line to line fault (AB), assuming the grounding resistance is zero, then the 
instantaneous boundary conditions will be: 


I; =0 l =- Ig and Va = Vg (4) 


Then, the boundary condition instantaneous will be: 
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Iz la 

Ip -1/3 v3 I, 

Io |= 0 (5) 
l, -la — 1/3 V3 h 


2.2.3. Line to Line to ground Fault in Phase AB to G 
A suppose line to line to ground fault (ABG), assuming the grounding resistance is zero, then the 
instantaneous boundary conditions will be: 


Ic = 0,1,= Ip and V= V= 0 (6) 


Then, the boundary condition instantaneous will be: 


lx 2/3 la— 1/3 I, 

Ip 1/3 V3 I, 

Io |= 1/31, +1/3 1, (7) 
ly —2/3 I,+1/3 I, + 1/3 V3 I, 


2.2.4. Three Phase Fault in Phase ABC 
A suppose for three phase fault (ABC), assuming the grounding resistance is zero, then the 
instantaneous boundary conditions will be: 


la + lp + le = 0 and V+ ý+ V=0 (8) 


Then, the boundary condition instantaneous will be: 


2 1 1 
Ty 3 la E 3 Ip J 3 I, 
Ip 243 i AN3 i 
3 3 
Ip] = r (9) 
L 2 1 1 a 1 1 
Y = latz l= 3V3 I + bt Nek 


Based on the above analysis, the characteristics of various faults based on Clarke’s transformation 
a- Modal, f- Modal and modal y can be proposed. 


2.3. Wavelet Transformation 

Wavelet transformation is the decomposition of a signal by a function @, çt) which is deleted and 
translated by the so-called mother wavelet. The function of the mother wavelet can be written as follows 
[22], [23]: 


ary) = 9 (—) (10) 
Where a is the dilation parameter ( a € Real ) and b is a translation parameter ( b € Real ), parameter a 
indicates the width of the wavelet curve, when the value of a wider magnified wavelet curve is diminished as 
the curve gets smaller, while the wavelet parameter curve b shows the localization of wavelet centered at 
t = b. The detection of the discrete wavelet transformed (DWT) fault is required so that the equation 
becomes [24], [25] 


Par(y=2!" p (2/(a—b)),j.keZ (11) 
Variables j and k are integers that scale the shifts of the mother wavelet function, to produce types of 


mother wavelet such as Syms and Haar wavelets. The width of a wavelet is shown by the scale a, and the 
position is indicated by the wavelet scale b. 


Algorithm for Fault Location and Classification on Parallel Transmission Line .... (Makmur Saini) 


702 o ISSN: 2088-8708 


Discrete Wavelet Transformations (DWT) are methods used to decompose the input signal, and the 
signal is analyzed by giving treatment to the wavelet coefficients. The decomposition process involves two 
filters, which are low-pass filter and a high-pass filter [26]. This is achieved by successive high pass and low 
pass filtering of the time domain signal and is defined by the following equation: 


ônign Ik] = Xn X[n]. g[2k - n] (12) 

Siow [k] = Èn X[n] . h[2k = n] (13) 
where Spign [k] = output high-pass filter and tow [k] = output low-pass filter. 

The signal is first passed through the high-pass filter and low pass filter, and then half of each output 


is taken as sampling, down through the sampling operation, the signal of desired frequency component can be 
obtained from the recurring decompositions as shown by Figure 1 [27]. 


gin] © 
ain © hn} t) 


Level 3 
hfn] Coefficient 
n 
"HE 
Level 2 
Coefficient 
"HO 
Level 1 
Coefficient 


Figure 1. The process decomposition of discrete wavelet transform 


This is called a decomposition first level process in the frequency range of 250-500 kHz. The output 
from the low-pass filter is used as the input of the next decomposition second level, with a frequency range of 
125-250 kHz. This process is repeated again until it reaches the third level of decomposition with a frequency 
of 62.5-125 kHz, and so forth, according to the level desired. In this study, it is chosen to reach level 4, with a 
frequency range of 31.25-62.5 kHz and a level 8 with frequency range of 7.8125-15.625 kHz [28]. 

The combination of the output from the high-pass filter and low-pass filter output is called 
Coefficient Wavelet transformed (CWT), which contains information on the results that have been 
compressed and transformed. In this study, the high pass filter and low-pass filter are coupled and becoming 
a Quadrature Mirror Filter (QMF), in which the couples meet the following [29] equation: 


hin] = (-1)" g[L+1-n] (14) 


where h[n] = high-pass filter, g[n] = low-pass filter, L= length of each filter. Successful to the down 
sampling operation that removes redundant information signal, wavelet transform has become one of the 
most reliable and accurate composition methods [30]. 


3. THE PROPOSED ALGORITHM 
In this study, the simulations were performed using PSCAD, and the simulation results were 
obtained from the fault current signal. 
The steps performed in this study were: 
a. Finding the input to the Clarke’s transformation and wavelet transform. The signal flow of PSCAD 
was then converted into m. files (*. M) and then converted into mat. Files (*mat) with a sampling 
rate of (10°) and frequency dependent of 0.5 Hz-1 MHz. 
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b. Determining the data stream interference, where the signal was transformed by using the Clarke’s 
transformation to convert the transient signals into the basic current mode signal in using 
Equation (1). 

c. Transforming the mode current signals again by using DWT and WTC, which were the generated 
coefficients, and then squared to be (WTC)? in order to obtain the maximum signal amplitude to 
determine the time of the interruption [31], [32]. 

d. Processing the ground mode and aerial mode and (WTC)’ using the Bewley Lattice diagram [32], 
[33] of the initial wave to determine the fault location and detection. 

e. Adding the magnitude of the current gamma to make the fault classification algorithm which then 

became modal J, that can be proposed. While the magnitude of the current gamma of each 

different types of fault can be seen in Table 1. 


I, =- 2/3 1,41/3(1+ v3) I, + 1/30-V3) L (15) 


f. Determining whether the internal fault in both Circuit 1 or Circuit 2 and external fault using the 
following equation: 


If a >1 or s >1 Fault Internal Circuit 1 (16) 
a2 B2 

If va |> 1 or e | >1 Fault Internal Circuit 2 (17) 
ad B1 

If faa |_ 4 or z =1 Fault External (18) 
Ta2 Tp2 


g. The protection technique should be able to classify the faulted phase for single-phase-to-ground 
faults. In the case of single-phase-to-ground faults, two of the modal components that include the 
faulted phase should have almost the same amplitude and the other modal component should be 
zero, as follows: 


lel — 1 | Le >AG fault (19) 

Y 

a -1 | Be >BG fault (20) 
Y 

ee 4 | Zé >CG fault (21) 
Y 


The algorithm will continue to determine the faulted phases involved in a multiple-phase fault. In 
the case of line to line faults, the criteria are as given in (22)-(24): 


lal [| = liyl| <o >AB fault (22) 
llil + |1|- Hal] < o >AC fault (23) 
llil- lhl] < o = BC fault (24) 


In the case of double line to ground faults, the criteria are as given in (25)-(27): 


(I+ — (lal + lib <6 >ABG fault (25) 
\(fal + Wold — Cie] + D| < 6 >ACG fault (26) 
(ly | + Hol) — Mal + [ie )| < 6 =BCG fault (27) 


h. Determining the fault Classification where the fault classification was divided into 2 categories 
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ground and unground. The current approach had zero given threshold or less equal to zero than the 

disruption of unground fault, otherwise if the current J) is greater than the specified threshold limit, 

it would be the ground fault. 

The unground fault was the line to line fault, while the threshold limit was given for termination 
criteria o = 0:02, while the ground was divided into 2, which were line to ground fault with the given 
threshold e = 0:03 and the line to line to ground fault with the given threshold 6 = 0:05 for the termination 
criteria [34], as shown in Figure 2 that illustrates the proposed fault-type classification algorithm that uses the 
modal components of the current signals. 


4. SIMULATION RESULTS AND DISCUSSIONS 

The system was connected with the sources at each end, as shown in Figure 2. This system was 
simulated using PSCAD/EMTD. For the case study, the simulation was modelled on a 230 kV double circuit 
transmission line, which was 200 km in length. 


Figure 2. One line diagram of the simulated transmission system 


Transmission data: 
Sequence Impedance ohm/km. 
Transmission LineZ; = Z2 = 0.03574 + j 0.5776 Zo = 0.36315 + j 1.32.647 


Fault Starting = 0.22 second Duration in fault = 
0.15 Second 
Fault resistance (Rp) = 0.001, 25, 50, 75 and 100 ohm 


Fault Inception Angle = 0, 15, 30, 45, 60, 90, 120 and 150 degree 
Source A and B Z; = Z, = Zo = 9.1859 + j 52.093 Ohm 

Type Conductor = Chukar, diameter = 1.602 inch, 

D, = 0.0524 ft = 0.0162763 m 


The results of the calculations took into account the position of the tower and distance between the 
conductors. The conductor types used for this simulation were obtained using the propagation 


velocity = aa = 299863.4379 km/second. 


4.1. Internal Fault 

Table 1 shows the results of various internal fault variations with fault resistance = 0.001 Ohm and 
the fault inception angle of 0 degrees at various distances of errors. The biggest fault currentoccurred on the 
line to ground fault (AG) at a distance of 25 km with a current of 2.8953 kA. The selected threshold line to 
ground for fault e = 0.03 from the simulation results showed that the fault threshold line to ground (AG) was 
e = 0.0237 which is smaller than the threshold (e) set. 

Also shown in Table 2 Ig,/Ig2 and Ig,/Ig2 on all types of faults that are greater than 1, which indicate 
that the fault was an internal fault in circuit 1, WTC of the aerial mode and ground mode are shown in 
Figure 3 The Clarke’s transformation and the mother wavelet Db8 obtained the biggest percentage error in 
fault location, about 0.4649 %, at the three phase disorder (AB) which had a distance of 75 km and the WTC 
of the aerial mode and ground mode are shown in Figure 4 while the smallest percentage error in fault 
location was about 0.0750 %, in line to line (AC) at a distance of 100 km. This indicates that the proposed 
algorithm for fault classification is accurate and precise. 
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Table 1. The obtained result for different faults using the Clarke’s Transformation and mother wavelet Db8 


Fault AG BG AB AC ABG ACG ABC 
Distance (km) 25 50 75 100 125 150 175 
Rs ( Ohm) 0.001 0.001 0.001 0.001 0.001 0.001 0.001 
Fault Inception Angle 0 0 0 0 0 0 0 
la (kA) 2.8953 0.5330 -2.3266 2.3941 -1.6213 1.5088 1.3156 
I, (kA) 0.5543 3.3782 2.6240 -0.4809 2.3185 -0.5580 1.7432 
I, (kA) 0.5590 -0.5358 0.4682 2.1064 -0.5121 -1.9317 -1.8289 
Tq, (KA) 1.9146 -0.8914 -2.3264 2.3938 -1.6364 1.6941 1.3141 
Ip, (kA) 0.4722 2.0446 1.6739 1.0421 1.4762 0.9851 1.8197 
la (kA) 1.9612 2.9020 4.0132 -1.3441 3.0248 1.2014 2.3410 
I, (kA) 1.0492 1.1841 0.00178 0.00198 0.5472 -0.4538 0.0321 
laz (KA) -0.5674 0.4403 0.6615 -0.5065 -0.6707 0.7786 -0.8424 
Igo (kA) -0.3011 -0.4040 0.4453 0.5377 0.6624 0.4012 0.9667 
Tor / 12 -3.3743 -2.0245 -3.5168 -4.7267 -2.4391 2.1758 -1.5599 
Ig, / Igo -1.5682 -5.0609 -3.5212 1.9381 -1.9417 2.1708 1.8823 
Calculate point of 25.765 49.103 75.930 99.850 123.995 150.98 174.235 
Fault (km) 
% Error 0.3825 -0.4994 0.4649 -0.0750 -0.5024 0.4994 0.3825 
x 10° Aerial Mode W” Ground mode 


T T T T 5 T T T T T T 


T %3 T T y T T T 
Hæ—— ta=0.22009 Second g-———Hta = 0.22009 second } 


o 0 
0.22 0.2201 0.2202 0.2203 0.2204 0.2205 0.2206 0.2207 0.2208 0.2209 0.221 0.22 0.2201 0.2202 0.2203 0.2204 0.2205 0.2206 0.2207 0.2208 0.2209 0.221 
Time (Second) Time (Second) 


(a) (b) 


Figure 3. Line to line fault (AG) circuit 1 from 25 from bus A, (a) Aerial mode for wavelet mother Db8, (b) 
For ground mode for wavelet mother Db8. 


x16" Aerial Mode x 107" Ground Mode 


822 0.2201 0.2202 0.2203 0.2204 0.2205 0.2206 0.2207 0.2208 0.2209 0.221 


i i i i i i i i i 
0.22 0.2201 0.2202 0.2203 0.2204 0.2205 0.2206 0.2207 0.2208 0.2209 0.221 
Time (Second) 


(a) (b) 


Figure 4. Line to ground fault (AB) circuit 1 from 75 from bus A, (a) Aerial mode for wavelet Mother Db8, 
(b) For ground mode for wavelet mother Db8 


4.2. Influence of Fault Resistance 

Table 2 shows the effects of variations in the resistance of 25 ohms and 50 ohms with fault inception 
angle at 15 and 45 degrees with varying distances. The simulation showed that when the resistance was 
increased, the fault current would drop, where the biggest fault current occurred in line to ground (BG) at a 
distance of 50 km with a fault resistance of 25 ohms. 
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Table 2 shows the effects of variations in the resistance of 25 ohms and 50 ohms with fault inception 
angle at 15 and 45 degrees with varying distances. The simulation showed that when the resistance was 
increased, the fault current would drop, where the biggest fault current occurred in line to ground (BG) at a 
distance of 50 km with a fault resistance of 25 ohms. The fault inception at the angle of 15 degrees was 
I, = 2.2829 kA and the threshold obtained at e = 0.00056 was smaller than the threshold set ¢ = 0.03. 
When using the mother wavelet Db8, the percentage error in the fault location was ranged between 
0.4-0.5%, as shown in Table 3. Iq4/Iq2 and Ig,/Ig2 varied from 2.5-8, indicating that the fault was internal to 
circuit 1. 


Table 2. Result for different resistance faults based on mother wavelet Db8 


Fault BG AB ACG ABC 

Distance (km) 50 75 125 150 

Fault Inception Angle 15 15 45 45 

Fault Resistance Rr (Ohm) 25 50 25 50 25 50 25 50 

la (kA) -0.2241 -0.2277 -.2.0771 -1.6669 1.4632 -1.1601 1.5970 1.4178 
I, (kA) -2.2829 1.7328 2.1951 1.7831 0.2377 0.2190 1.8558 -1.576 
I, (kA) 0.2450 -0.2375 0.1574 0.1574 -1.4154 -1.1325 -1.7960 1.5244 
Tq, (KA) -0.9056 -0.7183 -2.0771 -1.6669 1.2068 -1.0092 -1.5970 1.4178 
Ip, (kA) 1.2291 0.9303 1.3429 1.0968 0.9234 0.7546 1.8761 1.5944 
Ty (kA) 2.1335 -1.6427 3.4108 2.7635 -0.7128 0.6131 2.5946 2.2336 
I, (kA) 0.7745 0.5730 0 0 0.4156 -0.2974 0 0 

Tq (KA) 0.2046 -0.2253 0.4731 0.4566 0.4330 0.3544 0.6073 0.5599 
Tgp (KA) 0.2581 -0.2329 0.3001 0.2724 0.4203 0.3759 0.5707 0.5207 
la /le -4.4262 3.1882 4.3904 3.6507 2.7871 2.8476 2.6296 2.5322 
Ig, 1 Igo 8.2661 3.9944 11.3655 4.0264 2.1970 2.0074 3.2874 2.1773 
Calculate point of Fault (km) 49.760 49.010 75.93 76.048 124.100 123.99 150.83 150.98 
% Error -0.1199 -0.4994 0.4649 0.5024 -0.4499 -0.5024 0.4199 0.4949 


4.3. Influence of Fault Inception. 

As shown in Table 3, the simulations showed the effect on the variation of fault inception angle, 
ranging from 30 degrees to 150 degrees, with variations in fault resistance of 75 ohm and 100 ohm in the 
various types of fault and fault distance. Meanwhile, the threshold obtained in line to line ground disturbance 
(BCG) at a fault distance of 75 km, fault resistance 75 ohm and fault inception angle 75 degrees was 
ô = 0.0851; greater than the threshold set on the fault line to line to the ground of ò = 0.05, as it was faulted 
to the ground. The fault line to line (AC) at a fault distance of 125 km, fault resistance of 100 ohms and fault 
inception angle of 30 degrees obtained a threshold of o = 0.0017, smaller than the threshold set of o = 0.02. 
Table 3 also shows that if the fault inception angle was enlarged, then the fault current would increase, except 
for fault three-phase (ABC), which resulted with I41/I«2 and Ip,/Ig2 , between 1.2-3, indicating that the fault 
was internal fault circuit 1. Meanwhile, the percentage error in fault location was obtained around 0.2-0.5% 
when using the mother wavelet Db4. The fault classification algorithms show that the proposed algorithm is 
correct and responsive. 


Table 3. The obtained result for different Fault Inception Angle based on mother wavelet 


Fault AG AC BCG ABC 

Distance (km) 150 125 715 50 

Fault Resistance (Ohm) 75 100 75 100 

Fault Inception Angle 60 75 30 45 75 90 120 150 

Tq (kA) 0.8616 0.9626 -0.8963 -0.9292 0.5147 -0.6369 -2.2589 -2.0813 
l, (kA) 0.3496 0.4948 0.0639 0.1791 1.4682 1.5235 2.1835 -1.9497 
I (kA) 0.4069 0.5048 0.9368 1.0783 1.4639 1.5305 2.1482 2.0210 
Tay (KA) 0.6723 0.8005 0.8963 0.9292 0.7293 -0.8503 2.2589 2.0813 
Tg, (kA) 0.3219 0.4762 0.5763 0.7183 1.4764 1.5586 -2.1178 -1.9726 
Ty KA) 0.7216 0,7418 -0.3217 -0.2940 1.7983 1.9133 3.2379 2.8001 
I, (kA) 0.1957 0.1829 0 0 0.3223 0.3496 0 0 

Tq (KA) 0.4196 0.5677 0.2908 0.3547 0.4812 0.6391 0.9405 1.0646 
Igo (kA) 0.2260 0.3706 0.2374 0.3752 0.5994 0.7006 1.0632 1.1361° 
Tay / laz 1.6002 1.4100 3.0822 2.6197 1.5165 -1.3305 2.4018 1.9555 
Ip, /Ipe 1.4187 1.2849 2.4276 1.9144 2.4631 2.2247 2.2517 1.7363 
I, (kA) 0.1957 0.1829 0 0 0.3223 0.3496 0 0 
Calculate point of Fault (km) = 150..39 150.54 124.59 124.29 75.63 75.71 49.46 49.01 
% Error 0..1949 0.2699 -0.2204 -0.3525 0.3150 0.3525 -0.2699 -0.4946 
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4.4. Percentage Error in Fault Location for different Types of Mother Wavelet 

Table 4 shows the simulation error in determining the fault location when using various types of 
mother wavelet on the interference resistance of 100 ohms, fault inception angle of O degrees and a distance 
of interference, varied from 25 km-190 km. The fault location estimation error was formulated as follows: 


Fault Simulated Distance—Real fault Distance 


x 100% (24) 


Error % = 
s Total lenght of line 


From Table 4 it can be seen that the smallest average error was shown by the mother wavelet DB4, 
while mother wavelets Sym 4, Coil 4 and Db8 have a variation around 0.2-05% 


Table 4. Percentage error in fault location for different types of Mother Wavelet from Circuit 1, 
Rs = 100 Ohm and Fault Inception Angle (0) = 0 


Type Of Actual Db4 Coif4 Syms4 DB8 

Fault Point Calculate Error = Calculate Error = Calculate Error = Calculate Error = 
Fault point of (x-d)/L point of (x-d)/L point of (x-d)/L point of (x-d)/L 

(km) Fault (km) *100% Fault (km) *100% Fault (km) *100% Fault (km) *100% 

25 25.0151 0.0076 25.7650 0.3852 25.3150 0.1575 25.7650 0.3825 

50 49.7601 -0.1949 49.6110 -0.4949 49.0103 -0.4949 49.0103 -0.4949 

75 75.2550 0.1275 75.8549 0.4274 76.0080 0.5024 76.0048 0.5024 
LG 125 124.745 -0.1275 124.142 -0.4274 123.995 -0.5024 123.995 -.0.5024 
(AG) 150 150.239 0.1949 150.389 0.4949 150.989 0.4949 150.988 0.4949 
175 174.985 -0.0076 174.350 -0.3852 174.695 -0.1575 174.235 -0.3825 

190 189.981 0.0091 189.082 -0.4590 189.750 -0.1216 189.232 -0.3840 

25 24.5650 -0.2174 24.1153 -0.4423 25.4651 0.2325 25.9150 0.4575 

50 49.7601 -0.1199 49.1603 -0.4199 48.8603 -0.5695 48.9350 -0.5323 

75 75.4050 0.2204 75.8549 0.4274 76.0080 0.5024 76.0048 0.5024 
LL 125 124.590 -0.2204 124.142 -0.4274 123.995 -0.5024 123.995 -.0.5024 

( AB) 150 150.239 0.1199 150.841 0.4199 151.139 0.5698 151.065 0.5323 
175 175.435 0.2174 175.885 0.4423 174.535 -0.2325 174.085 -0.4575 

190 190.432 0.2159 190.232 -0.3840 189.680 -0.1591 189.150 -0.4215 

25 25.0151 0.0076 25.0151 0.0076 25.0151 0.0076 25.6150 0.3075 

50 49.7520 -0.1237 49.3120 -0.3449 49.1610 -0.4199 50.8849 -0.4425 

75 75.6290 0.3150 75.8849 0.4424 75.8540 0.4274 76.0048 0.5024 
LLG 125 124.371 -0.3150 124.161 -0.4224 124.145 -0.4274 123.995 -.0.5024 
(BCG) 150 150.241 0.1237 150.691 0.3449 150.839 0.4199 149.115 0.4425 
175 174.985 -0.0076 174.985 0.0076 174.985 -0.0076 174.385 -0.3075 

190 190.282 0.1409 189.682 -0.1591 189.982 -0.0091 189.240 -0.3765 

25 25.0151 0.0076 25.3150 0.1575 25.6150 0.3075 25.9890 0.4950 

50 49.4600 -0.2699 49.0103 -0.4949 48.8310 -0.5848 49.0103 -0.4949 

75 75.4050 0.2025 76.0048 0.5024 76.0080 0.5024 76.1090 0.5549 
LLL 125 124.595 -0.2025 123.995 -0.5024 123.995 -0.5024 123.901 -.0.5549 
(ABC) 150 150.540 0.2699 150.989 0.4949 151.169 0.5848 150.988 0.4949 
175 174.985 -0.0076 174.685 0.1575 174.395 -0.3075 174.020 -0.4950 

190 190.582 0.2909 189.532 -0.2340 189.530 -0.2340 189.090 -0.4515 


4.5. External Fault 
Table 5 shows that I«1 / Ia2 = 1 and Ig1/lg2 = 1 in various types of fault, with a fault resistance of 
100 and 0.001 Ohms, and fault inception angle of 0 and 60 degrees, which shows that the disturbance was an 
external fault. The largest fault current was found in the type of fault Bus B Line to line (AB) with fault 
resistance 0.0001 Ohm, fault inception angle at 60 degrees, and at 1.4344 kA. 


4.6. Determining Fault Location using Karenbauer’s Transformation 

Table 6 shows the percentage error in the calculation of fault location by using another method, 
which is Karenbaur’s Transformation, for comparison when using Clarke’s transformation (Table 1). It 
turned out that the results obtained by both methods in determining the fault location percentage error were 
similar. The only difference between Iq1/la2 and Ig,/Ig2 in, Table 1 was that the transformations achieved 
Tq / Tag = -3,3743 and Ig4/Ig2 = -1.5682 respectively, in the type of line to ground disturbance (AG) at fault 
distance 25 km, fault inception angle 0 degrees and fault resistance 0.001 ohm, whereas Table 6 shows type 
of fault at Igq /Igg = 3.6515 and Ig, / Igo = 2.4853. 
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Table 5. The obtained result for different external fault 


Fault Location Bus A Bus B 

Fault AG AB BCG ABC AG AB BCG ABC 
Fault Inception Angle 0 60 0 60 0 60 0 60 

Ry (Ohm) 0.001 100 0.001 100 

la (kA) -1.1360 1.2477 -0.3563 -0.6316 -0.7879 -1.4344 0.3333 -0.8333 
l, (kA) 0.3369 -1.5471 0.5576 -0.5190 0.3972 1.1107 0.3329 0.8671 
le (kA) 0.4560 0.3269 0.5570 0.6439 -0.3662 0.3271 0.3319 -0.8220 
Tay (KA) -0.9213 1.2477 -0.4162 0.6316 -0.6139 -1.4344 0.3228 -0.8333 
Tg, (kA) 0.3047 -1.0665 0.5654 -0.5342 0.3045 0.4572 0.3412 0.8856 
Ta (KA) -0.9213 1.2478 0.4166 0.6316 -0.6138 -1.4342 0.3231 0.8332 
Igo (KA) 0.3051 -1.0661 0.5658 0.5343 0.3048 0.4573 0.3415 0.8854 
Tay /la2 1 1 1 1 1 1 1 1 

Ip, 1 Ipe 1 1 1 1 1 1 1 1 

I, (kA) -0.2151 0 0.0624 0 0.2136 0 0.0631 0 

la (kA) -1.1360 1.2477 -0.3563 -0.6316 -0.7879 -1.4344 0.3333 -0.8333 


Table 6. The obtained result for different faults based on Karenbauer transformation and mother wavelet Db8& 


Fault AG BG AB AC ABG ACG ABC 
Distance (km) 25 50 75 100 125 150 175 

Ry (Ohm) 0.001 0.001 0.001 0.001 0.001 0.001 0.001 
Fault Inception Angle 0 0 0 0 0 0 0 

Ia (kA) 2.8953 0.5330 -2.3266 2.3941 -1.6213 1.5088 1.3156 
Ib (kA) 0.5543 3.3782 2.6240 -0.4809 2.3185 -0.5580 1.7432 
Ic (kA) 0.5590 -0.5358 0.4682 2.1064 -0.5121 -1.9317 -1.8289 
Tuy (KA) -1.0604 -1.0156 -1.6500 0.8938 -1.2098 0.6492 -0.9158 
Ig, (kA) 0.8711 -0.2779 -0.6763 1.5000 -0.4432 1.1054 0.9397 
la2 (kA) -0.2904 0.2989 -0.3918 -0.2586 0.2967 0.3248 0.5674 
Igo (KA) 0.3505 0.1752 -0.2992 -0.3757 0.3797 0.4649 0.5864 
Tay 1 Tae 3.6515 -3.3978 4.2133 -3.4563 -4.0776 1.9988 -1.6140 
Ip, 1 Ipe 2.4853 1.5862 2.2603 -3.9925 -1.1673 2.3778 1.6025 
Io (kA) 1.0422 1.1841 0.0179 0.0198 0.5472 -0.4538 0.0321 
Calculate point of 25.765 49.103 75.930 99.850 123.995 150.98 174.235 
Fault (km) 

% Error 0.3825 -0.4994 0.4649 -0.0750 -0.5024 0.4994 0.3825 


5. CONCLUSION 

The application of parallel transmission lines requires a more special consideration in comparison 
with the single transmission line, due to the effect of mutual coupling method. To overcome this problem, a 
new method was proposed, by using the Current alpha and beta (Current mode) from the Clarke’s 
transformation to convert the signal. Then, discrete wavelet transform (DWT) is used to obtain wavelet 
transform coefficients (WTC) ? to determine the current time when the fault amplitude values (WTC ? would 
reach a maximum point, to determine the fault location distance. This paper also proposed algorithm fault 
classification by using Clarke’s transformation. The simulation results showed that the results were accurate, 
which were also compared against the results obtained by using the Karenbauer’s transformation. 
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